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Abstract. The specific heat of the two-dimensional ±J Ising model has been 
investigated by the numerical transfer matrix method and Monte Carlo simulations 
from a new point of view. The region where a part of the specific heat takes the 
negative value has been investigated, which is characteristic of frustrated systems and 
reflects the non-trivial degeneracy of the ground state. The region mentioned above is 
found to be fairly large in the p — T plane (p is the concentration of the ferromagnetic 
bond and T is the temperature). Moreover, it includes the Nishimori line. Namely, 
it includes a part of the paramagnetic-ferromagnetic phase boundary, on which the 
specific heat cannot diverge. The present result indicates that the specific heat does 
not diverge at least on a part of the paramagnetic-ferromagnetic phase boundary above 
the multicritical point, which is in conflict with previous results. 



PACS numbers: 75.50.Lk,02.70.Lq,64.60.Cn,05.50.+q 
1. Introduction 

To elucidate the nature of critical phenomena of two-dimensional disordered Ising 
models has been a subject of a long-standing interest. Harris[l] concluded from a 
heuristic argument that the nature of critical phenomena of a disordered system becomes 
different from that of the corresponding pure system when the critical exponent of the 
specific heat, a, of the pure system is positive, while it remains the same when a < 0. 
Since the two-dimensional pure ferromagnetic Ising model is the marginal case, namely 
a = 0, many authors have investigated the properties of critical phenomena of the 
two-dimensional disordered Ising models [2- 14]. 

For the two-dimensional unfrustrated random Ising models, many authors 
concluded that the specific heat diverges double logarithmically at the paramagnetic- 
ferromagnetic phase boundary[2-9], though there are several results which insist that 
the specific heat remains finitefl, 10-12]. 

For the two-dimensional ± J Ising model which corresponds to the two-dimensional 
frustrated random Ising model, a few results exist. Inoue[13] investigated the system- 
size dependence of the height of the peak of the specific heat at p > 0.94, and 
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concluded that the logarithmic divergence of the specific heat at the paramagnetic- 
ferromagnetic phase boundary is most probable in the same way as the pure case. 
Reis et al[14] investigated the system-size dependence of the correlation length at 
p > 0.92, and concluded that logarithmic corrections do not play a role in contrast 
with the unfrustrated disordered systems. They concluded, however, that the specific 
heat diverges at most logarithmically. Namely, they could not completely exclude the 
possibility of the double logarithmic divergence of the specific heat. 

On the other hand, there exisits the Nishimori line, which is defined by the equation, 
exp(2J//c B T) = p/(l —p). (&b is the Boltzmann constant, which we put that A; B = 1 
from now on.) On the Nishimori line, several rigorous results have been derived [15]. 
Particularly, it has been proved that the specific heat remains finite. 

From the renormalization group approach, it is generally believed that, for the 
two-dimensional ±J Ising model, there exist two fixed points on the paramagnetic- 
ferromagnetid phase boundary above the multicritical point, namely the pure 
ferromagnetic fixed point and the multicritical fixed point[16-17], and there seems to 
be no random fixed point in contrast with the three-dimensional case [18]. (In this 
paper, we use the word "multicritical point" as the crossing point of the ferromagnetic- 
nonferromagnetic phase boundary and the Nishimori line, though there seems to be 
no spin glass phase in the two-dimensional case.) The pure ferromagnetic fixed point is 
stable, and the multicritical fixed point is unstable. Thus, the critical phenomena on the 
paramagnetic-ferromagnetic phase boundary above the multicritical point is governed 
by the pure ferromagnetic fixed point. Therefore, the critical exponent, a, should be 
zero on the whole paramagnetic-ferromagnetic phase boundary above the multicritical 
point. The fact mentioned above, however, does not make a help to determine the 
critical behaviour of the specific heat, since each of a logarithmic divergence, a double 
logarithmic divergence and a cusp like behaviour belongs to the case, a — 0. 

In this paper, we investigate the specific heat from a new point of view. We devide 
the specific heat of a system into two parts, C\ and C2. (For the detail definitions of C\ 
and C2, see section 2.) The value of C\ is easily found to be non- negative and remain 
finite at finite temperature. On the other hand, the value of C2 becomes negative at 
T = when the system is frustrated and the ground state has non-trivial degeneracy. 
Since we can consider that the negative value of C2 reflects the non-trivial ground state 
degeneracy, it is an interesting problem to make it clear up to what temperature the 
property persists. It is noted that in order that the specific heat may diverge, C2 should 
become infinite. 

Thus, we investigate the region in the p — T plane where C2 takes the negative value 
for the two-dimensional ±J Ising model. Our result shows that the region mentioned 
above is fairly large in the p — T plane. Moreover, it includes the Nishimori line. Namely, 
the region includes a part of the paramagnetic-ferromagnetic phase boundary above the 
multicritical point. In that region, the specific heat cannot diverge. 

The present result indicates that the specific heat does not diverge at least on a 
part of the paramagnetic-ferromagnetic phase boundary near and above the multicritical 



The specific heat of the two-dimensional ± J Ising model 



3 



point. The result is not directly in conflict with previous results[13,14] which insist that 
the specific heat diverges at most logarithmically at the paramagnetic-ferromagnetic 
phase boundary near the pure ferromagnetic case(p > 0.92), since the range of the 
concentration of the calcultions does not overlap. It is natural, however, to think that 
the nature of the phase transition does not change on the paramagnetic-ferromagnetic 
phase boundary above the multicritical point. Thus, we insist that the specific heat 
of the two-dimensional ±J Ising model does not diverge on the whole paramagnetic- 
ferromagnetic phase boundary above the multicritical point, though we cannot exclude 
the possibility that there exists a singular point on the paramagnetic-ferromagnetic 
phase boundary which separates the nature of the specific heat. 



2. The model and the method 



We consider the two-dimensional ± J Ising model on a square lattice with only nearest 
neighbor interactions. The Hamiltonian is written as follows: 

H = -^JijSiSj, (1) 

(ij) 

where Si = ±1, and the summation of (ij) runs over all the nearest neighbors. Each 
is determined according to the following probability distribution: 

P(J l3 ) = P 5(J tJ - J) + (1 - p)5(J ij + J). (2) 

In this paper, we put that J = 1. 

The specific heat per bond, C(p,T), is written as follows; 

C(P,T) = ^A[< E > T } P = ^^Et< -Wis, > T ] P , (3) 

where < • • • >t denotes the thermal average in a given bond configuration, {Jij}, at 
temperature, T. [• ■ -] p denotes the configurational average at the ferromagnetic bond 
concentration, p, and Nb is the number of bonds. 

Now, we divide the specific heat of the system, C(p, T) into two parts, Ci(p, T) and 
C 2 (P,T): 

C(p,T)=C7 1 (p,r) + C7 2 ( P ,r), (4) 

where 

1 d 

Cl(p,T) = J^Y,gf-i < -Ji3°i°i >{T}]p \{T}=T, (5) 
B (y) ij 



and 



1 d 

C 2 (P,T) = Tr^ ^ Wf^^ - J iJ a i a J >{T}]p \{T}=T ■ (6) 



Here, we introduce technically, Tjj, the local temperature of the bond, Jy. Namely, 
C\(p,T) is considered to be the configurational average of the change of the local 
energy, — < Jij^aj >{t}, when we infinitesimally increase the correspondeing local 
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temperature, T^. On the other hand, C 2 (p,T) is considered to be the configuratinal 
average of the the change of the local energy, — < JijOiOj >{t}, when we infinitesimally 
increase the temperarture, 7] m , which surrounds the local bond, J^. 
It is easily calculated that 



Therefore, at finite temperature, in order that the specific heat may diverge, C 2 (p,T) 
should become infinite. 

Now, we consider the case at zero temperature. At zero temperature, the specific 
heat, C(p, 0) becomes zero. When, there is no non-trivial degeneracy in the ground state 
of the system, each [< a^aj >\} p = 1. Namely, Ci(p, 0) becomes zero. Thus, C 2 (p, 0) 
also becomes zero. On the other hand, when the ground state has non-trivial degeneracy, 
some of [< OiOj >t\p become less than one. Namely, Ci(p, 0) becomes positive infinite. 
Thus, C 2 (p, 0) becomes negative infinite in this case. The negative value of C 2 (p,T) 
at finite temperature may be considered to be one of the influence of the non-trivial 
ground state degeneracy of the frustrated system. Thus, it is an intersting problem to 
make it clear up to what temperature the property persists. Therefore, in the following 
sections, we investigate the region in the p — T plane where C 2 (p, T) takes the negative 
value for the two-dimensional ±J Ising model. 

3. Results by the numerical transfer matrix method 

In this section, we investigate the region in the p — T plane where C 2 (p, T) takes the 
negative value by the numerical transfer matrix method. We have calculated for the 
lattice size, L = 4 — 16, at the bond concentration, p = 0.65 — 0.95. In the calculations, 
we take the free boundary condition, and take the configurational average for 10 4 — 10 5 
bond configurations. 

Figure 1 shows the temperature dependence of C(p,T), Ci(p,T) and C 2 (p, T) for 
L = 4 at p = 0.8. We can see that C 2 (p,T) takes the negative value in the low 
temperature region. We have estimated the temperature, T (L,p), where C 2 {p,T) takes 
the value zero for L = 4 — 16 at p = 0.65 — 0.95. Figure 2 shows the temperature 
dependence of C 2 (p ) T) near the zero point for L = 14 at p = 0.9. From the figure, we 
have estimated the temperture of the zero point, T (L = 14, p = 0.9) = 1.2764(14). The 
estimated values of T (L,p) for various L and p are shown in table 1. The accuracy of 
the values of T (L,p) becomes worse as the concentration, p, becomes small, since the 
gradient of C 2 (p,T) near the zero point becomes small. 

We can see that, at each concentration, p, the size dependence of the values of 
T (L,p) is very small, and the values are fairly large compared to the value of T N (p), 
the tempearature of the Nishimori line at p. 



(ij) 

Namely, C\ (p, T) has an upper bound; 




(7) 



Ci(p,T)<-L. 



(8) 
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Table 1. The estimated value of T (L,p). T^(p) is the temperature of the Nishimori 
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L 


p = 0.65 


p = 0.7 


p = 0.75 


p = 0.8 


p = 0.85 


p = 0.9 


p = 0.95 


4 


3.88(6) 


2.85(5) 


2.222(16) 


1.803(8) 


1.4968(28) 


1.2553(13) 


1.0303(20) 


6 


3.88(6) 


2.845(35) 


2.220(12) 


1.8045(45) 


1.5025(25) 


1.2672(12) 


1.0517(12) 


8 


3.88(12) 


2.835(65) 


2.234(22) 


1.8065(85) 


1.5020(40) 


1.2722(22) 


1.0650(20) 


10 


3.90(10) 


2.885(45) 


2.230(12) 


1.8035(75) 


1.5045(30) 


1.2744(26) 


1.0650(20) 


12 


3.88(4) 


2.83(5) 


2.220(14) 


1.8035(65) 


1.5038(25) 


1.2759(19) 


1.0823(18) 


14 


3.88(4) 


2.862(42) 


2.225(14) 


1.8055(35) 


1.5050(25) 


1.2764(14) 


1.0880(20) 


16 


3.89(7) 


2.845(30) 


2.225(15) 


1.803(7) 


1.5047(19) 


1.2768(18) 


1.0932(12) 


T N (P) 


3.2308 


2.3604 


1.8205 


1.4427 


1.1530 


0.9102 


0.6792 



Table 2. The estimated value of T (p) for various p. T^(p) is the temperature of the 
Nishimori line at p. 



p 


Up) 


Tn(p) 


0.65 


3.88(6) 


3.2308 


0.7 


2.845(65) 


2.3604 


0.75 


2.217(19) 


1.8205 


0.8 


1.805(9) 


1.4427 


0.85 


1.506(3) 


1.1530 


0.9 


1.278(3) 


0.9102 


0.95 


1.154(3) 


0.6792 



We have no definite principle to estimate the value of T (p) which is the extrapolated 
value of T (L,p) to L — > oo. No clear size dependence, however, can be seen at 
p = 0.65 — 0.8. Thus, we perform naive extrapolation to L — > oo in this region. It can be 
seen that the value of T (L, p) slightly increases as L increases at p = 0.85—0.95. We have 
found that the extrapolation by the iV _1 -law works fairly well at p — 0.85 and p = 0.9, 
where N is the total number of the spins of the system. Figure 3 shows a plot of T (L,p) 
versus 1/N at p = 0.9, where we have estimated that Tq(p = 0.9) = 1.278(3). We have 
also found that the extrapolation by the L _1 / 2 -law works fairly well at p — 0.95. The 
estimated values of T (p) for various p are shown in table 2. Strictly speaking, we cannot 
justify the above extrapolaton at each p. We can say, however, that the extrapolated 
values, T (p), might not change drastically even if we use other extrapolation laws, 
and there seems to be no possibility that the value of T (p) becomes smaller than the 
temperature of the Nishimori line at each p. 

Figure 4 shows the estimated values of Tq(j>) for various p in the p — T plane, which 
are denoted by black squares. The dashed line denotes the Nishimori line, and the open 
circles denote the paramagnetic-ferromagnetic phase boundary above the multicritical 
point [19]. We can see that the region where C^p, T) takes the negative value is fairly 
large. We can also see that the region mentioned above contains the Nishimori line. 
There are several numerical calculations about the paramagnetic-ferromagnetic phase 
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Table 3. The condition of the Monte Carlo simulation. 



p 


L 


M.C.S. 


N b 


0.88 


31 


400000 


240 




61 


800000 


96 




121 


800000 


120 


0.89 


31 


600000 


240 




61 


5000000 


120 




121 


16000000 


80 


0.9 


31 


800000 


320 




61 


1400000 


160 




91 


30000000 


64 


0.91 


31 


400000 


480 




61 


10000000 


96 




121 


20000000 


32 


0.92 


31 


600000 


720 




61 


1400000 


120 




121 


800000 


64 



boundary of the two-dimensional ±J Ising model [13, 19-22,27]. It is noted that all 
the calculations are almost consistent with each other and they show that, above 
the multicritical point, the concentration, p, of the phase boundary increases as the 
temperature increases. Thus, a part of the paramagnetic-ferromagnetic phase boundary 
is included in the region where C^p, T) takes the negative value. In that region, the 
specific heat cannot diverge. Namely, the specific heat does not diverge at least on a 
part of the paramgnetic-ferromagnetic phase boundary above the multicritical point. 

4. Monte Carlo simulations near the paramagnetic-ferromagnetic phase 
boundary 

In order to elucidate the region where C^p, T) takes the negative value near the 
paramagnetic-ferromagnetic phase boundary, we have performed the Monte Carlo 
simulations for more larger lattices, L = 31 — 121, at p = 0.88 — 0.92. In the 
calculations, we take the skew boundary condition in one direction and the periodic 
boundary condition in the other direction. 

The condition of the simulation is shown in table 3. M.C.S. denotes the Monte 
Carlo step of the simulation, which is chosen to be 20r, where the relaxation time, r, is 
evaluated by the statistical time-independent method[23], and Nb is the number of the 
bond configurations. 

The temperature dependence of C^Qo, T) near the zero point for L = 91 at p = 0.9 
is shown in figure 5, where we have estimated that T (L = 91, p = 0.9) = 1.2765(45). 
The estimated values of T (L,p) for various L and p are shown in table 4. 

Though T (L,p) takes similar value even if the lattice size, L, changes, we have 
extrapolated the values of T (p, L) to L — > oo by the iV _1 -law. We show the estimated 
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Table 4. The estimated values of T n (L,p) for various L and p. 

p L = 31 L = 64 L = 96 L = 121 

0.88 1.3485(55) 1.3535(45) 1.3555(35) 

0.89 1.3075(45) 1.303(5) 1.3055(55) 

0.9 1.286(4) 1.2775(45) 1.2765(45) 

0.91 1.269(4) 1.2665(45) 1.2685(45) 

0.92 1.249(4) 1.2485(45) 1.2485(35) 



Table 5. The values of T (p) for various p. T N (p) is the temperature of the Nishimori 
line at p. 

p T (p) T N (p) 

0.88 1.355(5) 1.0038 

0.89 1.3045(55) 0.9566 

0.9 1.275(5) 0.9102 

0.91 1.267(5) 0.8644 

0.92 1.2485(55) 0.8188 



values of T (p) for various p in table 5 with the temperature of the Nishimori line, T^(p). 

The results are also shown in figure 6. The black circles and squares denote To(p) 
evaluated by the Monte Carlo simulations and by the numerical transfer matrix method, 
respectively. The dashed line denotes the Nishimori line, and the open circles denote 
the paramagnetic- ferromagnetic phase boundary above the multicritical point [16]. It 
can be seen that both results by the numerical transfer matrix method and Monte 
Carlo simulations are consistent with each other. 

The crossing point of the paramagnetic-ferromagnetic phase boundary and the 
boundary of the region where C 2 (p, T) takes the negative value has been estimated 
to be 0.8985(15). Therefore, we conclude that the specific heat cannot diverge on the 
paramagnetic-ferromagnetic phase boundary at least for p c < p < 0.8985(15), where 
p c is the concentration of the multicritical point. There are many numerical estimates 
oi p c : 0.8905(5) [14],0.8872(8)[24],0.886(3) [25], 0.8906(26] and 0.8907(2) [27]. Recently, 
there is a conjecture about the exact value of p c , which insists that p c = 0.889972 [28]. 

5. Conclusions 

We have investigated the property of the specific heat of the two-dimensional ±J Ising 
model by the numerical transfer matrix method and Monte Carlo simulations from a 
new point of view. We have estimated the region where C^ip, T) takes the negative 
value in the p — T plane, which is one of the characteritic of the frustrated system, and 
reflects the non-trivial degeneracy of the ground state. By the numerical transfer matrix 
method for L = 4 — 16, and the Monte Carlo simulations for L = 31 — 121, the region 
have been found to be fairly large in the p — T plane, which includes the Nishimori-line. 
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Moreover, the region includes a part of the paramagnetic-ferromagnetic phase 
boundary, on which the specific heat cannot diverge. Thus, our results indicate that the 
specific heat cannot diverge on the paramagnetic-ferromagnetic phase boundary at least 
near and above the multicritical point. On the other hand, there are several literatures 
that, near the pure ferromagnetic point(j» > 0.92), the specific heat diverges at most 
logarithmically [13, 14]. Both results are not directly in conflict with each other, since the 
range of the concentration of the calculations does not overlap. It is natural, however, 
to think that the nature of the phase transition does not change on the paramagnetic- 
ferromagnetic phase boundary above the multicriticsl point. In the literatures which 
insist the divergence of the specific heat, the change of the peak height of the specific heat 
of various lattice size was investigated [13, 14]. It is a very subtle problem whether the 
specific heat diverges or remains finite, when the divergence is so weak as the logarithmic 
divergence. Our extrapolation, however, is straightforward and the values of Tq(p) might 
not change drastically even when there exist other system-size corrections, and there 
seems to be no possibility that the value of T (p) becomes smaller than the temperature 
of the Nishimori line at each p. Thus, we insist that the specific heat of the two- 
dimensional ± J Ising model does not diverge on the whole range of the paramagnetic- 
ferromagnetic phase boundary above the multicritical point, though we cannot exclude 
the possibility that a singular point exists on the paramagnetic-ferromagnetic phase 
boundary, which separates the nature of the specific heat. 
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Figure 1. The temperature dependence of C(p,T), Ci(p,T) and C2(p,T) for L = 4 
at p = 0.8. 
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Figure 2. The temperature dependence of C2(p,T) near the zero point for L = 14 at 
p = 0.9. 




Figure 3. A plot of T (L, P) versus 1/N at p = 0.9. 
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Figure 4. The estimated values of Tb(p) for various p in the p — T plane, which are 
denoted by black squares. The dashed line denotes the Nishimori line, and the open 
circles denote the paramagnetic-ferromagnetic phase boundary above the multicritical 
point. 




Figure 5. The temperature dependence of C2(p,T) near the zero point for L = 91 at 
p = 0.9. 
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Figure 6. The estimated values of T (p) for various p in the p — T plane. The 
black circles and squares denote the values of T (p) by the Monte Carlo simulations 
and by the numerical transfer matrix method, respectively. The dashed line denotes 
the Nishimori line, and the open circles denote the paramagnetic-ferromagnetic phase 
boundary above the multicritical point. 



